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ABSTRACT 

Radial stellar migration in galactic discs has received much attention in studies of galactic dynamics and chemical 
evolution, but remains a dynamical phenomenon that needs to be fully quantified. In this work, using a Tree-SPH 
simulation of an Sb-type disc galaxy, we quantify the effects of blurring (epicyclic excursions) and churning (change of 
guiding radius). We quantify migration (either blurring or churning) both in terms of flux (the number of migrators 
passing at a given radius), and by estimating the population of migrators at a given radius at the end of the simulation 
compared to non-migrators, but also by giving the distance over which the migration is effective at all radii. We confirm 
that the corotation of the bar is the main source of migrators by churning in a bar-dominated galaxy, its intensity 
being directly linked to the episode of a strong bar, in the first 1-3 Gyr of the simulation. We show that within the 
outer Lindblad resonance (OLR), migration is strongly dominated by churning, while blurring gains progressively more 
importance towards the outer disc and at later times. Most importantly, we show that the OLR limits the exchange of 
angular momentum, separating the disc in two distinct parts with minimal or null exchange, except in the transition 
zone, which is delimited by the position of the OLR at the epoch of the formation of the bar, and at the final epoch. We 
discuss the consequences of these findings for our understanding of the structure of the Milky Way disc. Because the 
Sun is situated slightly outside the OLR, we suggest that the solar vicinity may have experienced very limited churning 
from the inner disc. 

Key words. Galaxies: formation — Galaxies: evolution — Galaxies: spiral — Galaxies: structure — Galaxy: stellar 
content 


1. Introduction 


Radial stellar migration in the galactic discs has been at- 
tracting increasing attention, including some theoretical 


Sellwood fc Binney|2QQ2 Minchev fc Famaey|2QlQ Daniel 


Wyse||2014|) or numerical works ([Brnnetti et al. [pony 

Minchev et al.||2011[ |2012b[ |Di MattL et al.||2013| IKoskar] 


et al. 2013 Vera-Ciro et al. 2014), and studies based on 


observations of the stars of the Milky Way (e.g . Haywood 
Schonrich fc Binney||2009 Yu et al.]|2012). 


2008 


In galactic discs where non-axisymmetric potential per¬ 
turbations such as bars or sp iral arms occur, stars can gain 


or lose angular momentum (Lynden-Bell & Kalnajs 1972 
e.g.), leading to outward or inward migration, respectively. 
The result of this mechanism is that stars can be found 
at a galactocentric radius differing significantly from their 
birth radius. Another reason for apparent radial migration 
is simply the nature of orbits in axisymmetric or nearly 
axisymmetric potentials: stars oscillate radially around a 
guiding radius, the amplitude of the oscillations increas¬ 
i ng wi th the radial velocity dispersion. [Sellwood fc Binney| 
(2002) studied the impact of transient spiral arms on the 
change in angular momentum of stars in N-body simula¬ 
tions. They found that stars nearly at corotation with a 
spiral pattern experience the largest changes in angular mo¬ 
mentum, larger than what can be experienced by stars at 
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the Lindblad resonances of the spiral pattern or in the rest 
of the disc. This prominent effect of corotation was con¬ 


firmed, in particular by Roskar et al. (2012), in N-body sim¬ 
ulations exhibiting several spiral patterns amongst which 
the strongest one causes large changes in angular momen¬ 
tum for stars around its corotation radius. Other studies 


(Di Matteo et al.||2013 

Brnnetti et al. 

2011 Kubryk et al. 

2013, e.g.). In the case where both a 
are present, the bar-spiral resonance 

bar and spiral arms 
overlap studied by 

Minchev & Famaey 

(2010); Minchev et al. (|2011, 2012b) 


larger changes in angular momentum. In addition, while 
angular momentum can occur at corotation without any 
increase in radial kinetic energy, resonance overlap is ex¬ 
pected to increase radial energy: stars that were originally 
on nearly circular orbits can thus have significantly more 
eccentric orbits once they have experienced a change in an¬ 
gular momentum that is due to a resonance overlap. 


As it redistributes stars over the disc of a galaxy, ra¬ 
dial migration has been invoked as a possible explanation 
to a number of unanswered questions of galactic and extra- 
galactic stellar populati ons: the upturn of t he mean age in 


the outer parts of discs ([Roskar et al. ( 

12008) and e.g Bakos 

et al. (2008); Zheng et al. ([2014[)), t 

he formation of the 

Galactic thick disc ([Schonrich & Binney 

2009 Loebman 

et al. 

2011), the dispersion in metallicity at a given age 
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observed on stars in the solar vicinity ( |Haywood| 2008 
Schonrich fc Binney||2009| . 


Each of these problems requires a different level of mi¬ 
grations. For instance, U-shaped age profiles would require 
that stars from the inner disc of a galaxy are massively 
present in the outer regions. Similarly, shifting the mean 
solar-vicinity metallicity towards higher values, as has been 
suggested to explain the flat ness of the age-metall icity re¬ 
lation at the solar vicinity (Loebman et al. 2011), would 
require massive migration of stars to the solar radius. 
Explaining the range of metallicities in the solar vicinity 
is a more subtle problem because it may require only very 
little migration as it depends on the local (6-10 kpc) metal¬ 
licity gradient. If the gradient is steep, the presence of high- 
metallicity stars found at the solar radius may be explained 
by blurring alone, without any significant churning. That 
the gradient is steep is suggested by the data, which now 
heavily indicate that the Sun is at the interface between an 
inner and an outer disc t hat each have different chemical 
properties, as advocated in Haywood et al. (|2013|), and con¬ 
firmed on more extended data by Nidever et al. (2014). The 
recent results found for stars at the solar radius seem thus 
to dispute that strong migration by churning is necessary 
to explain the metallicity distribution. 

A common prediction of N-body models of galaxy evo¬ 
lution is the significant migration of inner dis c stars in the 


external parts of galaxy discs (see for example Roskar et al 
’2008); Minchev et al. (2011[ 2014)). In the case of the outer 


isc of the Milky Way, which contains stars with guiding 
radii greater than ^ 10 kpc - some of them with pericen- 
tres small enough to reach the solar vicinity- there is no 
evidence that this extreme migration has ever occurred at 
any time in the last 10 Gyr. At every age, stars of the 
outer disc are more metal poor than inner disc stars and 


more alpha-enhanced (Haywood et al. 
cates that the outer Galactic disc must 


tiave 


2013]), which indi- 
bllowed a dif- 


causes the increase of the azimuthal and radial frequen¬ 
cies, and subsequently an increase of the bar pattern speed. 

In addition, the characteristics of the thick disc of the 
Milky Way hardly require any significant migration, and 


is rich in gas and turbulent 

([Haywood et al.[[2013i Lehnert 

et al.||2014 Haywood et al.j 

[2015p. On the theoretical side. 


phenomenological radial migration is now implemented in 


several Galactic chemical evolution m odels (Schonrich & 
Binney 2009 Kubryk et al. 2014a|b ), there is still much 


ferent chemical evolutionary path (Snaith et al. 2014). Even 
without access to age information, large-scale spectroscopic 
surveys like APOGEE are confirming the substantial differ¬ 
ence between the inner and the outer Milky Way disc, as is 
shown by the properties of their stars in the [(a/Fe]-[Fe/H] 
plane or by the change of the metallicity dist ribution func¬ 
tion with distance from the Galaxy centre (| Anders et al.j 
|2014 Nidever et al.|20i4 ). It is striking in this context that 
regardless of the nature and formation of the outer disc of 
the Milky Way, stars from the inner thin disc have not been 
able to significantly pollute the outer parts of the Galactic 
disc, when models of radial migration would predict this as 
a common scenario. This also casts some doubts on the idea 
that migration due to the present bar/spiral is at the origin 
of the U-shaped age profiles observed in external galaxies, 
at least in bar-dominated systems where this inversion oc¬ 
curs outside of the outer Lindbland resonance (OLR) of the 
bar. 

However, it is possible that the first pattern that 
occurred in the galaxies had a lower speed than the 
current speed and that the OLR was at larger radii, which 
would have extended the churning action at larger radii. 
During the mass assembly of galaxies, the concentration 
of mass increases, and consequently the pattern speed 
increases, imp lying a smaller OLR radius with time. The 
simulations of Bournaud & Gombes (2002), which included 
gas accretion, showed such an increase in the bar pattern 
speed that is due to a central mass accumulation that 


uncertainty about the importance of the role of radial 
migration. For example, what fraction of stars are subject 
to migration and on which distances? Is the migration 
episodic or continuous? How are the different parts of 
the disc affected by migration? While tentative answers 
to these questions have been given in past works, from 
the simple mo del of recurring spi r al per turbations in an 
isolated disc in Sellwood fc Binney (|2QQ2|) to re c ent ni odels 
in a cosmological context by Minchev et al. (2014), the 
detailed quantitative effects of radial migration are still 
subject to debate. 


In this paper, we are interested in distinguishing be¬ 
tween the effects of ‘blurring’, or in other words, the radial 
migration that is due to epicyclic excursions around a fixed 
guiding radius, and ‘churning’, which is the radial migration 
that is due to a ch ange in this guiding radius ( according to 
the terminology of Schonrich & Binney (2009|). In partic¬ 
ular, we seek to quantify the fraction of the stars involved 
in the migration, the extent of the migration, the fraction 
of migrators at a given galactocentric radius, and the re¬ 
gions affected by the process. The study is performed using 
an N-body+SPH simulation of an Sb-typ e galaxy that was 
first presented in Halle & Gombes ( 2Q13| ). A few results on 
radial migrat ion for this simula t ions w ere presented in the 
Appendix of Di Matteo et al. (2014). In the present pa¬ 


per, we first briefly list the parameters of the simulation 
in Sect. 2. In Sect. 3 we study the density resonances, and 
in Sect. 4 we quantify the effects of blurring and churn¬ 
ing, their strength, and their extent. Section 5 discusses 
the behaviour of migrators near the outer Lindblad reso¬ 
nance, and in particular its role as a barrier for migrators. 
Finally, the kinematic characteristics of migrating stars are 
discussed in Sect. 6, and we conclude in Sect. 7. 


2. Numerical simulation 


We used one of the simulations presented in [Halle fc] 
Gombes (2013). We briefly summarize the main character- 
istics of the simulations. The y were perform ed with the 
N-body SPH code Gadget-2 (Springel 20051) and include 
stochastic star formation following a Schmidt law, kinetic 
feedback from core-collapse supernovae and some detailed 
cooling of the gas down to 100 K, with the possibility of 
including cooling by molecular hydrogen, whose local mass 
fraction is computed based on a semi-analytic recipe from 
Krumholz et al. (2008 2009); McKee & Krumholz (2010). 


The simulated Sb type galaxy includes 
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Table 1. Sb galaxy parameters 


Sb M,. 

Md 

Mf, 

Mg 

Th 

ad 

hd 

Tb 

ag 

hg 

[Mq] 

[Mo] 

[M®] 

[Mo] 

[kpc] 

[kpc] 

[kpc] 

[kpc] 

[kpc] 

[kpc] 

1.7 10^^ 

4.5 10^“ 

“TTTcP” 

0.9 10^'^ 

12 

5 

0.5 

1 

11.8 

0.2 


a stellar (s) and a gas (g) discs that both have 
Miyamoto-Nagai density profiles: 


Ps,g{R,z)— 

as^gR? + (^as,g + (ag^g + 


R? 


+ (as,s + ^Jz^ + hlg^ j [z^ + hig) 


( 1 ) 


a stellar bulge (6) and a dark matter halo {h) with 
Plummer density profiles: 


/ N 3 M 5 
phAn = - 7 —— 




( 2 ) 


The masses of the different components and density pro¬ 
files parameters are shown in Table There are initially 
400 000 particles of each component with particle masses of 
moM, and m^ for the dark matter, stars, and gas, listed 
in Table The softening length e used in the gravitational 
force computation and the smoothing length h used for hy¬ 
drodynamics are also listed in Table Particle velocities 
were a ssigned in the following way (see |Halle & Combes 
(2013) for details): For gaseous and stellar disc particles, 
we computed circular velocities from gravitational acceler¬ 
ations and applied an analytic asymmetric drift correction 
to obtain a more realistic profile. Radial velocity disper- 

^ ^ , 3.36GQS(r) ^ 

sions were derived as cr^(r) = ---, with nir) and 

K[r) 

S(r) the epicyclic frequency and the surface density at ra¬ 
dius r, respectively, and adopting a Toomre parameter Q 
equal to 1. The azimuthal velocity dispersions cre{r) were 
computed from the radial velocity dispersions cr^(r) using 
the epicycle approximation, and the vertical velocity disper¬ 
sion CFz{r) was derived by assuming isothermal equilibrium 
for the discs. For the spheroidal components, the veloc¬ 
ity dispersion is isotropic and was derived from the second 
moment of the Jeans equation. It can be noted here that 
our initial conditions for the evolution of the disc are not 
repr esentative of clumpy di scs observed at high-redshifts 
(e.g. Elmegreen et al. |2QQ7 ). We focused on the effect of 
a bar on the evolution of an already formed thin disc. In 
clumpy discs with a chaotic evolution, resonant patterns 
causing significant churning are difficult to form because of 
the higher velocity dispersion, and the radial migration is 
more likely to be some blurring. The uncertainty on the 
morphology and kinematics of galactic discs before bar for¬ 
mation is large. In addition, we did not take into account 
any thick-disc component that may have contributed sig¬ 
nificantly to the total mass of the disc at the ep och of bar 
formation in the Milky Way, for example (e.g. [Haywood 


et al.||2Q13 Di Matteo et al.|(2Q14 Snaith et al.|[2Q14| ). This 

kinematically hot component may have been less affected 


aywo 


Table 2. Resolution of the simulation 


e 

h mDM 

m* 

nig 

[pc] 

[Mo] 

[Mo] 

[Mo] 

100 

> 3.7 10® 

1.4 10® 

2.5 10"^ 


by ra dial migration (Loebman et al. 2011 


2011 ). 


Brunetti et al. 


The simulated time-span we used is 9 Gyr, which is sim¬ 
ilar to the estimated time-span during whic h the Milky Way 


has been evolving with no major mergers (Haywood et al. 
Hammer et alj20Q7 e.g.). Figureshows the surface 


2013 


density maps of ditterent components m the disc (gas, all 
stars, and new stars formed during the simulation). A bar is 
formed during the first Gyr. In this simulation, the feedback 
efficiency is 10% and there is no molecular hydrogen cool¬ 
ing. With these parameters, the surface density of gas and 
stars is relatively smooth because this feedback efficiency 
and the absence of molecular hydrogen cooling prevent gas 
and stellar clumps from forming. We note that compared t o 
many previous studies (for example Minchev et al.| (2011); 

dis 


Di Matteo et al. (2013|)), the stellar and gaseous discs are 
very extended (initial conditions were generated with an 
initial cut at R = 36 kpc for both gas and stars). This 
choice also allows studying regions outside the OLR at all 
times, which in this model is located between 11 and 23 kpc 
from the centre (see next sections). The maximum extent of 
discs is difficult to determine observationally, be it at high 
or low redshifts, which means that considering an extended 
disc is not necessari l y a sin gular case. The recent study by 
van Dokkum et al. (2014) showed, for example, that the 
disc of MlOl, previously strongly underestimated, extends 
to as much as 18 scale lengths in radius. Moreover, we here 
focused rather on the effect of radial migration with respect 
to the location of the main resonances than on the abso¬ 
lute sizes. Figure shows the azimuthally averaged surface 
density of the simulated disc components as a function of 
the galactocentric radius. The old stellar disc dominates the 
surface density at all radii throughout the simulated time- 
span. Its surface density profile is very stable, with only 
mild departures from the initial conditions, which is due to 
the presence of the bar and spiral arms formed during the 
simulation. 


3. Resonances in the disc 

The disc forms a central bar that persists during the whole 
simulated time-span. Transient spiral arms are also present. 
We determined the pattern speeds through a classic Fourier 
method: 

— At times spaced by a constant interval of 10 Myr, we 
performed a spatial Fourier transform of the surface 
density E(R, 6>) of the stars to obtain the dominant az- 
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Fig. 1. Surface density maps of the gas (left), stars formed during the simulation (middle), and all the stars (right). 
Face-on view boxes have a size of [80 kpc x 80 kpc] and edge-on views a size of [80 kpc x 20 kpc]. The colour scale, 
shown in the top right corner, is the same for all the views. 
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Fig. 3. Spectrograms of the m = 2 Fourier mode. Top left plot: integration on 9 Gyr from 0.5 Gyr to 9.5 Gyr. Other 

tv ( 

plots: integrations on 1 Gyr centred on the times specified in the plots. The black curves are fl{R) (solid), Q{R) --— 


(dot-dashed), ft{R) 


k{R) 


(dashed). The horizontal line represents the estimate of the bar pattern speed. The vertical 


lines represent the estimates of the bar ILR (blue), corotation (red), and OLR (green) radii. 



Fig. 2. Time evolution of the azimuthally averaged surface 
density of the old stellar disc without the bulge (solid), the 
young stellar disc (dashed), and the gas disc (dot-dashed), 
as a function of the distance from the galaxy centre. 


imuthal modes in the different radial bins. 


^{R) oc / : 

Jo 


^{R,0)e^^^d0. 


( 3 ) 


— We then performed a time Fourier analysis of the differ¬ 
ent modes in each radius bin to study their azimuthal 


speed. 


rtf 

TmiR,Uj)(X / 

Jti 


( 4 ) 


In the spatial Fourier transforms, the m = 2 mode usu¬ 
ally dominates, which corresponds to 7r-periodic patterns: 
the central bar or two-armed spirals. The contours of the 
obtained power in the R-Vl plane for this mode are shown 
in Fig. for a time Fourier-integration performed on the 
total time intervals of 9 Gyr (first panel) and of 1 Gyr cen- 

UJ 

tred on 1, 3, 5, 7, and 9 Gyr. Q = —, with m = 2 here, is 


m 

the pattern speed. The integration on 9 Gyr shows that the 
contours at low radii, which correspond to the bar, fill a Ct 
range from 10 to 30 km/s/kpc. The time interval of 
1 Gyr on which the time Fourier transform was then per¬ 
formed around specific times was chosen so as to optimise 
the determination of the pattern speeds: for a given time 
resolution, a too short integration time prevents correctly 
determining low pattern speeds and yields a poor resolu¬ 
tion in frequency, while a too large integration time leads 
to a determination of overly averaged pattern speeds if they 
change significantly during this time. Figure shows that 
the bar slows down with time. The bar transfers angular 
momentum to the rest of the disc and to the stellar bulge 
and d ark matter halo, as already d iscussed, among oth ¬ 
ers, by Debattista fc Sellwood (|1998|); Athanassonla ( 2QQ2|); 


Martinez- TO nesta et al.| (|2QQ6|); |Saha et al.| ( 2QT^ ; |Sah^ 

& Naab| (|2013P; |Di Matteo et al. (|2014p. The angular mo¬ 
mentum transfer between the different components of the 
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galaxy is detailed in Fig. At the beginning of the sim¬ 
ulation, the angular momentum is contained in the disc 
because the DM halo and stellar bulge have no initial rota¬ 
tion. In 9 Gyr, the disc loses 8% of its angular momentum, 
which is transferred mainly to the DM halo, while a slight 
fraction is given to the bulge. 


The contours of Fig. show the main peaks of power in 
the Q-R plane. The bar (power peak at low radii) is clearly 
the strongest 7r-periodic surface density perturbation with a 
rigid-body rotation in the integration periods of 1 Gyr, but 
spiral arms (peaks at larger radii) are visible as well. For ex¬ 
ample, the panel corresponding to t = 3 Gyr indicates spi¬ 
ral patterns rotating at 18 km/s/kpc and 7 km/s/kpc. 
These resonant modes with pattern speeds such that some 
of the resonance radii coincide might be sustained by non¬ 
linear mode coupling ( Tagger et aL|1987 e.g.). Some higher 
m-periodic features are also present (m = 3,4...) with var¬ 
ious pattern speeds, but they are weaker than the bar in 
the period from 1 to 9 Gyr. 


/ d^0 _3 # 

^ RdR 


. The corotation resonance (GR) 


z=0 


occurs for D = Dp, the inner Lindblad resonance (ILR) for 
D — Dp = — , and the OLR for D — Dp = — —. From a 


number of time Fourier transforms on 1 Gyr around times 
separated by 100 Myr, we estimated the values of the ILR, 
CR, and OLR radii of the bar as a function of time (see 
Fig.i. The bar slow-down causes the bar corotation radius 
to increase and allows the bar to become more elongated. A 
measure of the bar strength using the same spatial Fourier 
decomposition of coefficients at each time is shown in 
Fig-U The bar strength increases until t = 3 Gyr, when it 
drops before increasing again as a res ult of the angular mo - 
mentum exchange with the DM halo ( | Athanassoula|[2QQ2 ). 
The drop in strength corresponds to the start of the buck- 
ling of the bar. This vertical instability, leading to a thick 
X-shaped bar, can be seen in the edge-on views of the new 
stars and the whole stellar component of Fig. 




35 
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E 


15 


10 
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Fig. 5. Time evolution of the ILR, CR, and OLR radii of 
the bar (left y-axis) and of the bar speed D (right y-axis). 

The various resonances are expected to locally affect the 
angular momentum of the stellar disc and thus to generate 
some radial migration. We study next the radial migration 
that occurs between 1 and 9 Gyr of evolution (after the 
initial building of the bar). 

4. Blurring and churning 


Fig. 4. Top: Time evolution of the vertical component of 
the angular momentum Lz for different galaxy compo¬ 
nents. Bottom: Time variation of the fraction of the total 
angular momentum contained in each galaxy component 
-Lz Lz t=o 

Lz tot t=o 


By determining the angular speeds corresponding to 
the main contours, we obtained the radii at which stars 
on circular orbits are in resonance with the main pat¬ 
terns. We computed the angular speed D from the poten¬ 
tial by D^ = —and the epicyclic frequency k by 
R dRz=o 


In this section, we quantify different types of radial stellar 
migration. We examine 


— migration in terms of galactocentric radius (hereafter 
simply radius), that is, the difference of radius between 
two different times. Blurring and churning can both 
cause a change in the galactocentric radius of a star, 
thus this definition captures both types of migrators 
and can be considered as the overall migration experi¬ 
enced by disc stars. The comparison between radii at 
two different times during disc ev olution is wide l y used 
in the literature, see for example Roskar et al. ( 2QQ8|); 
Brunetti et al. 


Matteo et al.| (|2Q 13|); [Knbryk et al.| (|2Q13|); 
Minchev et al. ~pQ12b 2Q14p. 
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— determination of the relative minima and maxima of the 
oscillatory evolution of the radius. 

— use of a linear fit between the relative minima on the 
one hand and the relative maxima on the other to obtain 
a local minimum radius Rminit) and a local maximum 
radius i?max(^) at time t. 

— definition of the guiding radius at t by the average of 

these two radii: {R){t) = ^ -^max(^) ^ 

Figure shows examples of this determination of the 
guiding radius for two random stellar particles. The posi¬ 
tions of the particles are centred on the centre of mass of 
the whole galaxy at each time-step for this analysis. This 
method directly provides the guiding radius and also the 
amplitude of the radial oscillations at any time t. 


Fig. 6. Time evolution of an estimate of the bar strength 
from the m = 2 coefficient of the Fourier decomposition 
of the stellar surface density. The bar strengthens between 
0 and 1 Gyr before weakening when a buckling instability 
occurs around 3 Gyr, and it then strengthens again. 


— migration in terms of guiding radius, or churning, that 
is, the difference of the guiding radius betwe en two dif¬ 
ferent times. This follows the nomenclature of lSchonrichl 


& Binney (2009). Stars that increase the amplitude of 
their radial oscillations over time but do not change 
their guiding radius, that is, stars that only experience 
blurring, are not considered as migrators according to 
this definition. 


We recall that our model is not intended to reproduce ei¬ 
ther the characteristics of the Milky Way disc - as an ex¬ 
ample, the pattern speed of our simulated bar is lower than 
the pattern speed measured for the Galaxy (Gerhard 2011) 
and as a consequence, the main resonances are located at 
much larger distances from the centre than those measured 
for the Milky Way - or of any other specific galaxy. As a 
result, the exact values of spatial or kinematic variations 
experienced by stars in the model are not directly applica¬ 
ble to any galaxy. However, the effects described below can 
be considered as representative of those that any typical 
bar-dominated galaxy probably experiences. 


4.1. Determination of guiding radii and amplitude of radial 
oscillations 

A stellar radius oscillates around a guiding radius. In the 
case of low eccentricities, it is possible to determine this 
guiding radius by finding the radius Rc at which a star has 
the same vertical component of the angular momentum and 
a circular trajectory: = Rc'^drciRc)^ where Ucirc is the 

circular velocity obtained from the potential. The ampli¬ 
tude of the radial oscillations around the guiding radius is 

expected to scale as where ctr is the radial velocity 

dispersion, and ctr and the epicyclic frequency k are both 
functions of the radius R. 

The output snapshots of our simulation are separated 
by 10 Myr, which is small enough to allow computing the 
guiding radius at any time t by a simple method: 



Fig. 7. Determination of the guiding radius for two random 
stellar particles. Blue lines: galactocentric radii R{t). Red 
lines: Rmin(^)- Green lines: Rmax(^)- Black lines: guiding 

radii {R){t) = + 


4.2. Overall migration compared with churning alone 

Our aim is to compare the overall population of migrating 
stars, by blurring and churning, with that of migrators by 
churning alone. For this, we first examined the distributions 
of the variations of radius (/guiding radius) as a function 
of the initial radius (/guiding radius), for different time- 
intervals between 1 and 9 Gyr of evolution. The top set 
of plots of Fig. shows the difference of radii of stellar 
particles at final time R and initial time R as a function of 
the radius at R, while the bottom set shows the difference of 
guiding radii as a function of the guiding radius at R , where 
R and R are given a range of values. Thus, in the top set of 
plots the whole population of migrators is shown, whilst in 
the bottom plots only migrators by churning are selected. 
The considered stellar components are the old stellar disc 
and the new stars that are formed before R. The RMS of 
the variations are indicated for each time interval. 

In each of the panels of Fig. the shaded areas show 
the span of the values of the ILR ^lue), CR (re^ and OLR 
(green) radii of the bar determined from Fig. at R and 
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Fig. 8. Distributions of the variations of galactocentric radii (top half) and guiding radii (bottom half). The colour of 
each bin signifies the stellar mass in the bin and the colour scale shown at the bottom of the figure applies everywhere. 
The shaded areas represent the bar ILR (blue), CR (red), and OLR (green) radii variation in the time-span of each plot. 
The red vertical lines are the average of the bar CR radius during the time-span of a plot. The diagonal red lines help 
identify the migrators that cross the bar CR radius from lower radii (dot-shaded part of the schematic plot) or from 
larger radii (grid-shaded part). The diagonal green lines similarly help identify migration with respect to the OLR radius 
at the end of the time interval. 
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tf. The vertical solid line is the average of the corotation 
radius and the diagonal line equation is ^ = Rcr — Ri for 
the radii plots, y = Rqr — {Ri) for the guiding radii plots. 
The regions between the diagonal and vertical lines allow 
estimating the migrators that cross the corotation radius 
by blurring+churning or by churning alone: 

— In the region corresponding to the dot-filled region of 
the schematic diagram, the stars in the radii plots have 
Ri < Rcr and Rf > Rcr, and those on the guiding 
radii plots have (Ri) < Rcr and (i?f) > Rcr- 

— In the region corresponding to the grid-filled region of 
the schematic diagram, the stars in the radii plots have 
Ri > Rcr and Rf < Rcr, and those on the guiding 
radii plots have {Ri) > Rcr and (i?f) < Rcr- 

Analogously, in each panel, the vertical dashed line shows 
the position of the OLR at t = tf, and the dashed 
diagonal line equation is ^ = Rqlr — Ri for the radii 
plots, y = Rolr — (^i) for the guiding radii plots. The 
regions between the diagonal and vertical dashed lines 
thus allow estimating the migrators that cross the OLR, 
at its final position in the time interval considered, by 
blurring+churning or by churning alone. We use these 
criteria in Sects. 1431 and [5l 

The results presented in Fig. can be summarised as 
follows: 


1 . 


As a general behaviour, the distribution for the whole 
population and for migrators by churning alone shows 
some features that appear as diagonal structures. These 
features are seeded by several resonances that are due to 
the presence of the bar and of the spiral arms. The most 
prominent diagonal features occur around the corota¬ 
tion of the bar, the main source of migrat o rs in the 


2 . 


stellar disc, as in Minchev fc Famaeyl ( |2Q1Q[ ); [Minchev 
(2011); Brnnetti et al. " pQlip 


et aL 


Quantifying radial migration by means of the instan¬ 
taneous difference between radii at two different times 
leads to overestimating churning both in terms of max¬ 
imum extent and RMS of the variations and in terms of 
fraction of migrators (see further discussion in Fig. 9). 
As an example, the RMS values for the 1 to 3 Gyr time- 
interval is 2.5 kpc for the variations of radius, while it is 
only 1.8 kpc for the variations of guiding radius (we note 
that the absolute values are specific to our model and 
may not be directly applicable to any specific galaxy, 
MW included). These differences are due to the radial 
excursions that can cause a stellar particle to have a 
radius between {R){t) — A{t) and {R){t) + A{t)^ where 
A{t) = Rmax{t) — {R){t) is the semi-amplitude of the 
radial oscillations at time t. The radius R{to) at time to 
can be expressed as 


R{to) = {R){to) + r{to), (5) 

where r(io) is in the interval [—^ 4 (^ 0 ), ^(to)]- The vari- 
ation in radius R between time U and tf is thus 


R(tf) - R{ti) = (R)(tf) - (R)(ti) + r(tf) - r(ti), (6) 

and r(tf) — r(ti) is in the interval [—(A(tf) + 
A(ti)), A(tf) + A{ti)]. 


Table 3. RMS values of the variations of radius induced by 
blurring alone, at different times. See Fig. for a compari¬ 
son with the variations obtained for the whole population of 
migrators and for those migrated by churning alone. Note 
that these values are obtained with an isolated galactic disc 
with specific size parameters (possibly more extended than 
the Milky Way disc) and are not applicable to any specific 
galaxy. 


II 

CO 

o 

tf = 5Gyr 

if = 7Gyr 

if = 9Gyr 

ti — IGyr 1.73 kpc 

1.73 kpc 

1.68 kpc 

1.95 kpc 

ti = 3Gyr 

1.84 kpc 

1.83 kpc 

2.04 kpc 

ti = 5Gyr 


1.90 kpc 

2.05 kpc 

ti = 7Gyr 



2.13 kpc 


3. The migration in terms of change of radius (=blur- 
ring+churning) keeps approximately the same distri¬ 
bution and RMS value on time-intervals of the same 
length, although the amplitude and RMS of the migra¬ 
tion in terms of change of guiding radius (=churning 
alone) are smaller at late times. For example, for time- 
intervals of 2 Gyr, the RMS value of the radius varia¬ 
tion only varies from 2.2 to 2.5 kpc, while the amplitude 
and RMS value of the guiding radius variation decreases 
from 1.8 kpc in a time interval of 2 Gyr from t = 1 to 
3 Gyr to 1.1 kpc between t = 5 and 7 Gyr or t = 7 and 
9 Gyr. This difference is due to the increase, with time, 
of the amplitude of radial oscillations, that is, blurring, 
as quantified in Table where the RMS of the distri- 
butiorQis given for the same time intervals as shown 
in Fig. In particular, from the comparison of the val¬ 
ues given in Table with those of churning alone given 
in the bottom panels of Fig. it can be noted that in 
the early phase of disc evolution, the spatial variations 
induced by churning overwhelm those due to blurring. 
This trend in reversed already after 3 Gyr of disc evolu¬ 
tion, when the extent of the radial variations by blurring 
dominates the variations induced by churning alone. 

4. The disc experiences its most intense phase of migra¬ 
tion by churning in the very early phase of its evolu¬ 
tion, when the stellar bar forms from an axisymmetric 
potential and it is thin and strong. After the buckling in¬ 
stability and the formation of the boxy/peanut-shaped 
bulge at t~3Gyr, the amplitude of the variation of the 
guiding radii of stars migrated by churning alone dimin¬ 
ishes. The increase in the RMS variation of the guid¬ 
ing radii as the time-span increases (cf., for example, 
the RMS variations of guiding radii in the time inter¬ 
val 1-3 Gyr versus 1-9 Gyr) does not necessarily imply 
that the number of migrators by churning increases with 
time, but only that migrators by churning have time to 
reach larger distances on longer timescales. We elabo¬ 
rate on this point below. 

To emphasize the difference between migration simply 
defined as a change in the instantaneous radius and mi¬ 
gration due to a change of guiding radius, the decrease of 
churning with time, and to quantify the fraction of stars 
involved in blurring and/or churning in this simulation, we 

^ The RMS values were evaluated as the quadratic difference 
between the RMS value of the overall population of migrators 
and that migrated by churning alone. 
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X[kpc] X[kpc] X[kpc] 


Fig. 9. Migration between 1 and 3 Gyr. First line: Mass 
fraction of stars that migrate by more than x kpc inwards 
or outwards in terms of radius (solid line) or guiding radius 
(dashed line). Second line: Mass fraction of stars that mi¬ 
grate outwards (right half of the plot) or inwards (left part) 
by more than x kpc in terms of radius (solid) or guiding 
radius (dashed). Third and fourth lines: Same study, but 
only for stars with a radius (solid lines) or a guiding radius 
(dashed lines) around Ri. 

represent in the top panels of Fig. the mass fraction of 
the stars that migrate by more than x kpc in terms of ra¬ 
dius or guiding radius from 1 to 3 Gyr. Figures [To] and [TT1 
then show the same fraction, but for stars migrating in the 
time interval between 5 and 7 Gyr and over the whole in¬ 
terval 1-9 Gyr, respectively. Figure for example shows 
that between 1 and 3 Gyr, 40% of the ensemble of stars 



Fig. 10. Migration between 5 and 7 Gyr. First line: Mass 
fraction of stars that migrate by more than x kpc inwards 
or outwards in terms of radius (solid line) or guiding radius 
(dashed line). Second line: Mass fraction of stars that mi¬ 
grate outwards (right half of the plot) or inwards (left part) 
by more than x kpc in terms of radius (solid) or guiding 
radius (dashed). The yellow lines in these two plots are the 
corresponding migration from 1 to 3 Gyr (identical to the 
curves of Fig. |^. Third and fourth lines: Same study, but 
only for stars with a radius (solid lines) or a guiding radius 
(dashed lines) around Ri. 

(with mass-weighting accounting for the different masses of 
stellar particles) change their instantaneous radius by more 
than 2 kpc, while this fraction reduces to only 20% for 
a change of the guiding radius of the same amplitude. The 
plotted fraction decreases more rapidly with x for the vari- 
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Fig. 11. Migration between 1 and 9 Gyr. First line: Mass 
fraction of stars that migrate by more than x kpc inwards 
or outwards in terms of radius (solid line) or guiding radius 
(dashed line). Second line: Mass fraction of stars that mi¬ 
grate outwards (right half of the plot) or inwards (left part) 
by more than x kpc in terms of radius (solid) or guiding 
radius (dashed). The yellow lines in these two plots are the 
corresponding migration from 1 to 3 Gyr (identical to the 
curves of Fig. |^. Third and fourth lines: Same study, but 
only for stars with a radius (solid lines) or a guiding radius 
(dashed lines) around Ri. 


ation of guiding radius and it extends to a lower maximum 
value. The plots in the second row show the detail of the 
inward and outward migration. In the negative x left part 
of the plot, the curves represent the fraction of stars mi¬ 
grating inwards by more than \x\ kpc, while in the right 
positive X part, the curves represent the fraction of stars 
migrating outwards by more than x kpc. The curves are 
almost symmetric with respect to the x = 0 vertical axis. 
The migration has a lower amplitude in terms of guiding ra¬ 


dius than in terms of radius for both inwards and outwards 
migration. The lower panels show the same study for stars 
in specific initial radius or guiding radius bins. More specif¬ 
ically, in these last two rows we quantify the stellar mass 
migrating of more than AR (showing both the absolute and 
the real value of the variation) from some initial radii i?i, 
normalising this stellar mass to the mass initially contained 
at = i?i. In particular, we chose to analyse the fraction of 
migrating stars with respect to three different initial radii, 
around Ri = 5 kpc, Ri = 10 kpc and Ri = 20 kpc, cor¬ 
responding to radii between the ILR and the GR, the GR 
and the OLR, and outside the OLR, respectively. From in¬ 
spection of these plots, we can deduce the following: 


1 . Over the time period 1-3 Gyr, the whole population of 
migrators moves limitedly with respect to its initial ra¬ 
dius, independently of the location of this initial radius 
with respect to the bar’s resonances. As an example, the 
fraction of stars that has a variation of its radius in ab¬ 
solute values greater than | AR| = 4 kpc is about 10% at 
Ri = b kpc, and 20% at Ri = 10 kpc and Ri = 20 kpc. 
These fractions decrease when considering stars that ex¬ 
perience churning, the fraction of stars that change their 
guiding radii by more than 4 kpc is smaller than 10% 
for all values of the initial radius inspected. 

2. If one considers the sign of the displacement (bottom 
panels in Fig. - with positive AR corresponding to 
outward migration and negative AR to inward migra¬ 
tion - , one sees that migration is, in general, not sym¬ 
metric. For radii between the ILR and the GR, there is 
an excess of stars migrating outwards, AR = —Ri be¬ 
ing the largest possible extent of inwards migration for 
migrators originating at R = Ri. Between the GR and 
the OLR, there is an excess of inward migrators, both 
in terms of blurring-fchurning and in terms of churning 
alone. We note in particular that outward migrators by 
churning do not cross the OLR at t=3 Gyr, the fraction 
of stars with AR > 4 kpc - this value is the smallest dis¬ 
placement required to cross the OLR during this time 
interval - is indeed null. A limited number of stars that 
originated between the GR and the OLR can cross the 
OLR by blurring, but this number constitutes a very 
limited fraction of the stars born inside the OLR (fewer 
than 5% for Ri = 5 kpc). We discuss this important 
finding in more detail in the next section. Finally, most 
of the stars that initially were located at Ri = 20 kpc, 
thus outside the OLR between 1 and 3 Gyr, do not reach 
the OLR, with migration mostly spatially redistributing 
stars that are ab initio in the outer disc. Moreover, there 
is a slight excess of inward migrators by blurring in this 
region. 


The results found in Fig.j^are confirmed when analysing 
migration over the same time duration (2 Gyr), but for a 


different time interval (5-7 Gyr). Figure 10 shows the same 
study as Fig. for the time period from 5 to 7 Gyr, with 
the bottom panels detailing the migration for initial radii 
located again between the ILR and the GR (Ri =7 kpc), 
between the GR and the OLR (Ri =15 kpc), and beyond 
the OLR (Ri =25 kpc). Because the bar slows down with 
respect to earlier times, the initial radii Ri are now more 
external than those chosen for the analysis shown in Fig.|^ 
The striking difference between migration between 5 and 
7 Gyr, compared to the time interval 1-3 Gyr (Fig. [9), is 
the even smaller change in guiding radius experience by 
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the stars at these late times: there are no migrators by 
churning with \AR\ > 4 kpc, most of the stars (80% for 
the region between the CR and the OLR, almost 100% for 
stars originating inside the CR or outside the OLR) expe¬ 
rience a change in their guiding radii smaller than 2 kpc. 
The strength of churning, in terms of spatial extent covered 
by the process, thus decreases. We also confirm what we al¬ 
ready noted in Figs. an di that is, migrators by churning 
originated in the inner disc, that is, inside the OLR, do not 
cross this resonance: the maximum displacement of outward 
migrators by churning whose Ri = lb kpc at = 5 Gyr is 
about 3.5 kpc, which is not enough to cross the OLR, whose 
location is at Rqlr = 20 kpc at time t = 7 Gyr. 


It is not straightforward to identify the reason of the 
decrease of the importance of churning with time that is 
observed in our modelled galaxy. It may be due to a con¬ 
comitance of factors that may be difficult to separate. This 
decrease with time can be explained by the evolution of the 
bar, and of its strength, the influence of other resonances 
in the disc, and the kinetic state of the disc. As an exam¬ 
ple, Fig.j^shows that after an abrupt decline, the bar takes 
time to regain strength in the final several Gyr. It is then 
not as concentrated in the galactic plane (it has a peanut 
shape), and the tangential force it exerts on stars in the 
plane is thus weaker than at earlier times. In the period 
from 1 Gyr to about 3 Gyr, there is also a phenomenon 
of r esonance overlap that can i ncrease the radial migration 
(see Minchev & Famaey (2010). T^ disc also gradually be¬ 
comes hotter (as we show in Sect.[^ see especially Figs. 16 


and 
perturbations. 


21 ), and stars thus become less responsive to potential 


Finally, when analysing migration over the whole time 
interval 1-9 Gyr, as done in Fig. [m one can see that the 
maximum displacement of migrators has increased. This is 
because stars have had time to cross a larger portion of the 
disc. For example, 5% of migrators by churning migrated 
by more than |AR| = 6 kpc, whilst this fraction was null 
between I and 3 Gyr. But when examining the details of mi¬ 
gration, it is still valid that the amount of the displacement 
depends both on its sign (inward versus outward migra¬ 
tors) and on the location of the initial radius Ri. We note 
in particular that the largest displacement of outwards mi¬ 
grators occurs in the region between the ILR and the CR 
{Ri = 6 kpc and Ri = 9 kpc, respectively), but that even 
over such a long time interval, these outwards migrators do 
not cross the final OLR radius. Furthermore, blurring sig¬ 
nificantly contributes to shaping the outer disc - compare 
for example the fraction of outwards migrators in the outer 
disc {Ri = 30 kpc) by churning+blurring versus churning 
alone. 


4.3. Local flux of migrators 

We now attempt to determine which migration effect dom¬ 
inates at some radius (migration only due to epicyclic ex¬ 
cursions away from a guiding radius or a change of guiding 
radius), before computing the fraction of different types of 
migrators at a final radius. 

A quantification of a ‘migration fiux’ as a function of the 
radius can be obtained by computing the mass fraction of 
the disc that crosses a given radius. We represent in the top 
plots of Fig. as a function of the galactocentric radius R 


— the whole migration flux, that is, the mass fraction of 
the stellar disc that crosses R, the crossing of R being 
quantified in terms of radius, that is: Ri < R, Rf > 
R from time U to time R for outward migration (solid 
lines) and Ri > R, Rf < R for inward migration (dashed 
lines) 

— the migration flux by churning alone , that is, the mass 
fraction of the stellar disc that crosses R in terms of 
radius and also of guiding radius: Ri < R, Rf > R and 
(Ri) < R, (Rf) > R for outward migration (solid lines) 
and Ri > R^ R{ < R and (Ri) > R, (Rf) < R for inward 
migration (dashed lines). 

Both fluxes in the top panels of Fig. [^are normalised to 
the total stellar mass in the disc. 

The whole outward migration flux at the corotation ra¬ 
dius is thus the portion of stars that belongs to the dot-filled 
area of the schematic diagram of Fig. in the Rf — Ri vs Ri 
plots, and the whole outward migration flux at a radius R 
different from the corotation radius is simply obtained by 
shifting the threshold radius and the diagonal line accord¬ 
ingly. Similarly, the whole inward migration corresponds to 
the grid-filled area of the schematic diagram. The migra¬ 
tion flux by churning alone is obtained by selecting stars 
that are both in the shaded zones of the radii plots and of 
the guiding radii plots (top and bottom panels in Fig. 
respectively). The migration flux by churning alone is thus 
always inferior to the whole migration flux because it is 
a fraction of it. The inner kpc inward migrators (both in 
terms of change in their instantaneous radius or in their 
guiding radius) are stars captured by the bar. Some other 
much smaller peaks are observed at large radii, correspond¬ 
ing to spiral patterns. However, in agreement with previous 
results (see Introduction), the most significant part of the 
migration flux occurs near the corotation of the bar because 
the bar is the strongest potential perturbation in our sim¬ 
ulation and its corotation radius is located in a disc region 
with a high surface density. In discussing these plots, we 
wish to emphasize two points: 

— The finding that the corotation is the locus of the 
strongest flux of migrators by churning can be partic¬ 
ularly appreciated considering the time intervals of dif¬ 
ferent lengths, for example comparing the time interval 
1-3 Gyr with those at 1-5 Gyr, 1-7 Gyr, and 1-9 Gyr. 
The corotation spans an increasingly larger radial ex¬ 
tent as the duration of the time interval grows. As a 
result, the flux through the corotation region, which is 
mostly a thin spike around the GR for short time inter¬ 
vals (cf. 1-3 Gyr), transforms into a sort of large plateau 
for longer time intervals, whose extent corresponds to 
the spatial extent spanned by the corotation during the 
corresponding time. 

— The dominant role of churning near corotation can also 
be appreciated in the bottom panels of Fig.[^ where we 
show the fraction of stars migrated by churning with re¬ 
spect to the whole sample of migrators crossing a given 
radius R . While more than 60% of migrators crossing 
the corotation region are migrating by churning, this 
fraction significantly decreases for other regions of the 
disc. In particular, it can be appreciated that the frac¬ 
tion of migrators by churning decreases for later times 
outside corotation. This confirms the growing impor¬ 
tance of radial heating (i.e. blurring) as time increases 
(see also discussion in the previous section). 
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Fig. 12. Top half: Mass fraction of the stellar disc crossing a radius R in terms of radius (apparent migration) or both 
radius and guiding radius (true migration). Bottom half: Ratios of the true to apparent migrators. The shaded areas 
represent the bar ILR (blue), CR (red), and OLR (green) radii variation in the time-span of each plot. 
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Fig. 13. Fraction of migrators as a function of radius (top half) or guiding radius (bottom half) at the end of the time- 
span specified in each plot. The shaded areas represent the bar ILR (blue), CR (red), and OLR (green) radii variation in 
the time-span of each plot. Note that these values are obtained with an isolated galactic disc with specific size parameters 
(possibly more extended than the Milky Way disc) and are not applicable to any specific galaxy. 
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4.4. Sources and sinks: How many migrators are there at a 
given radius? And how far from it did they originate? 

In the previous sections, we have quantified the main 
sources of migration in the disc and the flux of migrat¬ 
ing stars as a function of distance from the galaxy centre. 
We now discuss where migrators are redistributed in the 
disc, and in particular, how many migrators can be ex¬ 
pected at different radii throughout the disc. To this aim, 
we represent in Fig. the fraction of stars of a radius bin 
around Rf at R that have migrated by more than n kpc 
(n =2, 4, and 6) in radius since t\ (top half of the fig¬ 
ure) or whose guiding radius has changed by more than 
n kpc since ti (bottom half). These fractions peak on both 
sides of the CR radius of the bar, consistently with the 
migrating fluxes of Fig. which peak at the bar CR ra¬ 
dius: radii Rf > Rcr receive stars migrating outwards from 
Ri < Rcr 7 while radii Rf < Rcr receive stars migrating 
inwards from Ri > Rcr- However, as already discussed pre¬ 
viously (see also next section), inward and outward migra¬ 
tors cannot cross the whole disc, and in particular, outward 
migrators originated around the CR cannot reach the outer 
{R > Rolr) regions. 

At first glance, the fraction of the whole population of 
migrators (churning-{-blurring) seems very high at certain 
radii (see top panels in Fig. f^: as an example, in the time 
interval 1-9 Gyr, 40% of the stars in the CR region are 
stars that have reached this region by migrating by more 
than 2 kpc from some inner radii by blurring or churning; 
in the OLR region between 40% and 50% of stars are out¬ 
ward 2 kpc-migrators; outside the OLR, the contribution of 
these stars to the local population increases nearly mono- 
tonically up to 80%. This high fraction is due to the ex¬ 
ponentially declining stellar surface density that causes in¬ 
wards migrators to easily constitute a large part of the stars 
at an external final radius. The population of more extreme 
migrators - those that have migrated by more than 6 kpc 
outwards by blurring or churning - are still important con¬ 
tributors to the local (= at a given Rf) stellar population: 
the fraction has a local maximum in the OLR region, at 
radii < Rolr(^ = R), where it peaks at ^ 25%, it decreases 
to about 10% at Rf = Rolr, to finally increase again nearly 
monotonically in the outer disc. However, when we examine 
the contribution of migrators by churning alone at a given 
radius (bottom panels of Fig. 13), we can notice some dif¬ 
ferent trends and absolute values. In particular, extreme 
(AR > 6 kpc) outward migrators by churning never exceed 
20 % of the local stellar population, with this value reached 
in the OLR region at radii below the final OLR radius. After 
this maximum, the fraction diminishes to a value of about 
5% at the final OLR and stays constant in the outer disc. 
For inward migrators, the highest contribution is reached 
at radii just beyond the ILR region, with the contribution 
of long-distance migrators by churning being very limited. 

To quantify how migration (by churning and/or blur¬ 
ring) can affect a given radius, the distributions of initial 
and final radii for stars at differe nt radial bins in the disc 
are often shown (see for example, ( Minchev et al.|2QT4l). To 
allow for a comparison of our results with |Minchev et rL 
(2014), we show in Figs. 14 and 15 t he results of an analysi s 
similar to that in Figs. 3 and 4 of Minchev et al. (2014). 
To keep some consistency with the rest of our analysis that 
studies migration in a fixed amount of time, we used as 
birth radii the radii of all stars at 1 Gyr. The top row of 



Fig. 15. RMS values of galactocentric and guiding radii at 
1 Gyr of stars with galactocentric radii in 2 kpc wide bins 
around R at 9 Gyr. The vertical dashed lines, as in Fig.pTl 
are the GR (red) and OLR (green) radii at 9 Gyr. The 
vertical solid green line is the OLR radius at 1 Gyr. 


Fig. shows the distribution of guiding radii at 9 Gyr 
of stars whose radii are in different bins (shown in shaded 
colours). The histograms, especially at large radii, tend to 
peak at smaller radii than the average radii of the shaded 
area, exhibiting the well-known asymmetric drift effect. The 
second row shows the distribution of guiding radii at 1 Gyr 
of the same stars (with guiding radii in the shaded bins at 
9 Gyr). Gomparison of the top and middle row shows that 
the distributions of guiding radii change between the initial 
and final times. The amplitude of the distribution s tend to 
decrease at all radial bins (this is also observed by [Minchev 


et al. ( 2Q14| )). The initial distributions of (Ri), skewed to¬ 
wards low radii, translate into distributions of (R/) that are 
much more symmetric. The differences between the distri¬ 
butions of initial and final guiding radii is evidence that 
migration by churning has occurred, and in particular, the 
increase of guiding radii with time indicates a predomi¬ 
nance of outward migration. The bottom line shows the 
galactocentric radii of the same stars at 1 Gyr. The stars 
clearly have migrated significantly. The distributions are 
very similar to that of the middle line, indicating that near 
the beginning of the simulation, at t =1 Gyr, most of the 
stars are in circular orbits. We compare the distributions of 
initial radii and guiding radii more quantitatively in Fig.[T5| 
by representing the RMS value of the initial galactocentric 
radii or guiding radii of stars in bins of final radius (we 
use ov erlappi ng bins of 2 kpc widt h as in [Minchev et al.j 
( 2Q14[ )). As in Minchev et al. (2014), we see these RMS val¬ 
ues, estimates of the widths of the histograms of Fig. 
are generally slightly higher for the initial galactocentric 
radii. The implications of Figs. p!4[andpR|on the mixing of 
the disc are discussed in more detail in the next section. 


5. OLR: a barrier for migrators 

In the previous section, we have reported that migrators 
by churning do not cross the OLR. In this section we 
develop this point further. 

As previously described, the bottom panels of Fig. [^ show 
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Fig. 14. Distributions of final guiding radii (top), initial guiding radii (middle), and initial galactocentric radii (bottom) 
of stars whose final galactocentric radii belong to bins of 2 kpc width. Each curve corresponds to the closest same-colour 
shaded bin of final galactocentric radii. The time interval is from 1 Gyr to 9 Gyr. The vertical dashed lines are the GR 
(red) and OLR (green) radii at 9 Gyr. The vertical solid green line is the OLR radius at 1 Gyr. In the two bottom lines, 
the histograms corresponding to stars in radius bins outside the OLR radius at 9 Gyr have been plotted with thicker 
lines (see next Section). The stellar masses M in bins of each histogram are divided by the largest mass M^ax contained 
in a bin of the histogram. 


the variation AR = (Rf) — (Ri) of guiding radii as a func¬ 
tion of their initial guiding radius Ri in the time interval 
[ti, tf]. The dashed vertical line indicates the position of the 
OLR at tf, that is, at the end of the time interval under 
consideration, and the diagonal dashed line the equation 
AR = RoLR(^f) ~ (^i)- is evident from these figures 
that at all times, stars with (Ri) < RoLR(^i) do not cross 
the position of the OLR at the final time tf, these stars are 
always below the line AR = RoLR(^f) ~ (^i)- Migrators in 
the OLR region in the time interval under consideration, 
that is, with RoLR(^i) < Ri < Rolr(R), can cross the 
OLR, but these are stars that at least at t = ti had guiding 
radii all beyond the initial OLR position RoLR(^i)- In 
other words, the only migrators by churning in a barred 
galaxy at any given time that can be found beyond the 
OLR are stars migrating from a region between the OLR 
at the initial time of bar formation and the current OLR 
position. The larger the variations of the pattern speed 
of the bar over time, and consequently, the larger the 


variations ARqlr = -Rolr(R) ~ -RoLR(^i) d fde OLR 
position, the larger will be the region of the inner disc 
where migrators crossing the OLR can originate from. In 
no case, however, can stars born at corotation cross the 
final OLR. Analogously, if we were to draw a diagonal line 
of equation AR = RoLR(^i) ~ (Ri) in the bottom panels of 
Fig-i this line separating stars that can cross the initial 
position of the OLR, we would find that no stars with 
(Ri) ^ Rolr(R) would migrate inward crossing the radius 
R = RoLR(ti). This suggests that the OLR region, that 
is, the region between the initial (i.e. at the epoch of bar 
formation) and final (i.e. current) position of the OLR 
is a transition region for a disc galaxy, the only region 
where migrators by churning can be exchanged between 
the inner and the outer disc. But neither stars with 
(Ri) < RoLR(^i) can cross the OLR, nor, in the opposite 
sense, stars with (Ri) > Rolr(R) can penetrate below the 
radius Ri = RoLR(^i)- Note that some more contamination 
of the regions around the OLR can be caused by stars 
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migrating by blurring (top panels of Fig. [^, even if, in 
general, these polluters originate in regions outside those 
spanned by the CR. The middle row of Fig. also shows 
the very limited mixing between the regions inside the 
initial OLR (at t =1 Gyr) and the regions outside the 
final OLR (at t =9 Gyr). There are no stars that migrated 
by churning from inside the initial OLR and, vice versa, 
inside the initial OLR, there are no stars coming from 
outside the final OLR. The lack of contamination of the 
disc region outside the final OLR radius by the regions 
inside the initial OLR radius can be seen from the extent 
of the distributions of the initial guiding radii of stars 
that are beyond the final OLR (plotted as thicker lines). 
For example, none of the stars ending up at the final 
time between 24 and 26 kpc (that is, just outside the 
OLR at the final time t =9 Gyr) come from radii below 
Rolr(^ =1 Gyr)=ll kpc, whereas stars belonging to the 
bins centred between 5 and 21 kpc, inside Rolr(^ =9 Gyr), 
all have distributions of {Ri) that can reach values as 
low as 3 kpc. Figure shows that RMS values of initial 
galactocentric or guiding radii of stars in a final galacto- 
centric radius bin are monotonically increasing from the 
innermost regions up to a radius inside the final OLR. 
After this maximum, the RMS values tend to decrease in 
the outer disc, indicating again that there is no accumula¬ 
tion of migrators from the inner disc outside the final OLR. 


To our knowledge, this result has not been pointed 
out before, and it recalls the similar barriers encountered 


by th e gas in a barred galaxy (Schwarz 1981 Gombes 


1988): similarly to stars, gas in the CR-OLR region can 


also gain angular momentum to reach the OLR position, 
at most. Similarly to stars (Fig. bottom panels), gas 
accumulates in the OLR region, Rut differently from the 
stellar component, gas is then able to dissipate part of its 
energy, shocks and forms rings of concentrated material 
at the OLR. Similarly, gas outside the OLR is not able to 
penetrate this resonance and access the inner region of the 
disc, until the barrier is removed (because of a si gnificant 
change in the bar properties, see for example |Gombes| 

(MR). 


The results we presented are based on the analysis of 
a single simulation, and one may have concerns about the 
generality of the findings discussed. It is difficult to compare 
with other N-body works presented in the literature, either 
because they usually lack information about the position of 
the OLR at initial times and only give its average position 
or position at the final time analysed or because the OLR 
is too close to the initial boundary of the simulat ed stellar 
discs to make robust predictions (see for example Minchev 


et al. (2012b); Di Matteo et al. (2013)). However, an inter¬ 
esting comparison can be made with test particle simula¬ 
tions: in these experiments, the pattern speed of the main 
asymmetry is kept fixed, and as a consequence, the posi¬ 
tion of the OLR does not change with time; the role played 
by the OLR in th ese cases can be discussed. A comparison 
with the work byjMinchev & Famaey (2010), for example, 
shows that when only a bar pattern is present in the disc 
(their Fig 2, left panels), the angular momentum exchange 
induced by the bar diminishes closer to the OLR and stops 
very clearly outside it, with a null angular momentum vari¬ 
ation at radii outside the OLR. Independently of the bar 
strength, their experiments confirm the robustness of our 


results: the angular momentum exchange is only allowed 
inside the OLR, and, in particular, the maximal variation 
of angular momentum produced at the GR (where most of 
the migrators by churning are generated) is not sufficient to 
allow those stars to cross the outer resonance. In the right 
panels of the same figure, a similar analysis is made for a 
single spiral pattern. Unfortunately, the OLR is outside the 
x-range shown in these plots, and as a consequence, is not 
possible to discuss where the exchange of ang ular momen¬ 
tum stops. T his can be seen in the work by |Sellwood~~fc 
Binney (2002), however. Their constrained model, contain¬ 


ing a unique spiral pattern, shows that the OLR limits the 
exchange of angular momentum, with any variation stimu¬ 
lated at the OLR being extremely local and concerning a 
very limited number of stars. Of course, the response of a 
disc to asymmetric perturbations changes when many per¬ 
turbations of similar amplitudes are sup erimposed, as is the 


case o f several tr ansient spiral pattern s ( Sellwood fc Binney 
(2002), but also|Roskar et al. (2008) anT Vera-Giro et al. 


(2014)) or of the overlap of a bar and a spira l pattern of 
similar strengt hs (see, for example. Fig. 3 of Minchev & 
Famaey|(|2010)). But these cases are different from the case 
discusseci here, where a unique main pattern is present (the 
bar), with possibly an overlap with weak spiral arms, es¬ 
pecially in the first 2 Gyr of evolution of the system. In 
this case, most of the angular momentum exchange occurs 
near the corotation of the bar, far from the OLR, and the 
contamination of the outer disc {R > Rolr(R) with stars 
migrated from the inner disc {R < Rolr(R) is very weak 
or null. In this case of a growing bar, the stopping of an¬ 
gular momentum exchange at the OLR of the bar is more¬ 
over consistent with the mechanism of bar growth: bars 
strengthen and grow by receiving angular momentum from 
stars near their ILR and by giving angular momen tum to 
stars near their OLR (Lynden-Bell & Kalnaj s|1972| ). There 
is thus some mixing close to the OLR, but beyond the OLR 
region there are no more resonant stars that may exchange 
angular momentum with the bar. 

Why is this result important? For the Milky Way, at 
least for two reasons: 


1 . It may explain why the inner and the outer discs of the 
Galaxy have been able to maintain two different stellar 
populations over a time interval of ^10 Gyr, which prob¬ 
ably corresponds to the whole interval of secular evolu¬ 
tion ex perienced by the Galax y. We recall here the re¬ 
sults of [Haywood et ah] (2013), who extensively showed 
and discussed that the chemical properties of outer-disc 
stars {{R) > 10 kpc) are substantially different from 
those of stars of the inner dis c of similar ages. T hese 
results have been con firmed by Anders et al. (2014) and 
Nidever et al. (2014), who also showed that the chem¬ 
ical properties of the outer disc are signif icantly differ¬ 


ent fro m those of the inner disc (see also Bensby et al. 
(2012), for a high-resolution spectroscopy study of the 


inner and outer disc). The current position of the OLR 
in the Milky Way is estimated to be slightly inside the 


2000 

Famaey et al. 

2005 

Minchev et al. 


2007| ). This may therefore ex plain why the Sun app ears 


to be in a transition region ( Haywood et al.|[2013 ) and 


why the outer disc does not show signatures of any sig¬ 
nificant pollution by stars originating in the inner disc 
{R< 6-7 kpc). 
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2 . Depending on the exact location of the Sun with respect 
to the OLR, the solar vicinity may have been strongly 
polluted by stars migrating by churning from the inner 
disc (see Fig.[^ bottom panels, for stars inside the final 
OLR position^ or not at all. If the Su n is outside th e 
OLR position, as has been deduced by Dehnen (2000); 


Famaey et al. (2005); Minchev et al. (2007) and the bar 


is the main source of asymmetric perturbations in the 
Milky Way, then the impact of churning at the solar 
vicinity may have been significantly overestimated (see 
Fig. bottom panels, for stars beyond the OLR re¬ 
gion). In particular, stars that migrated by churning at 
the solar vicinity may well originate in regions much 
closer to the solar radius than previously thought, and 
in this case, invoking substantial migration from the 
corotation and the inner disc is not possible. 


This result may also lead to doubts about the interpreta¬ 
tion of U-turn in age profiles or inversion in colour-profiles 
found in the outer disc of external galaxies. When these in¬ 
versions occur outside the OLR position for bar-dominated 
galaxies, it is difficult to explain them in terms of strong 
migration from the inner disc, since, according to our re¬ 
sults, no strong migration from the corotation is expected 
to contaminate the outer disc, unless another pattern with 
lower speed existed before the present one, with an OLR 
extending to outer radii. 


6. Migration and cooling/heating 

6.1. Vertical heating/cooling 


The effect of radial migration on the vertic al structure of 
discs has been studied in a number of work s (|Minchev et al. 


2Q12al [Solway et "ar]|2Q12| [Roskar et al.||2Q13[ |Vera-Ciro 


et al. |2Q14[ e.g.). The assumption that outward migration 
would help create a thick disc because of the higher verti¬ 
cal velocity dispersion of the outward migrators originating 
from the inner hotter disc has been debated because of the 
vertical cooling outward migrators should undergo if the 


vertical action of a stellar orbit is conserved ([Minchev et al. 


2012a Solway et al.||2Q12 Roskar et al.||2Q13 )^ 


We observe that the vertical velocity dispersion of the 
stellar disc slightly increases with time at all radii. In Fig.p^ 
we plot the vertical velocity distribution of the old stellar 
disc, the stellar disc of stars formed during the simulation, 
and the gas disc. The velocity dispersion of the new stars, 
formed from the colder gas disc, is generally lower than the 
velocity dispersion of the old stellar disc. The mass of new 
stars is only a small fraction of the total stellar mass (see 
Fig. however, so star formation has no net cooling effect 
on the total stellar disc. 

We would like to know if the stars that are going to 
migrate are a special sub-population in terms of kinetics 
of stars with the same initial guiding radius, whether the 
migration affects their kinetic state, and if the migrators 
aff ect the population of st ars at their final guiding radius. 

Vera-Ciro et al. ( |2Q14 ) found a provenance bias of mi- 
grators in terms of kinetic state. They found that radial 
migration driven by spiral arms (seeded in their simula¬ 
tions by massive perturbers) preferentially affect stars with 
a lower velocity dispersion than the average velocity disper¬ 
sion at the initial radius. We find a similar trend, as shown 
in Fig. where the ratio of the vertical velocity disper¬ 
sion of stars in bins in (Ri) and (Rf) — {Ri) (we only used 



Fig. 16. Time evolution of the vertical velocity dispersion 
of the old stellar disc without the bulge (solid), the young 
stellar disc (dashed), and the gas disc (dotted-dashed). 


bins that contained at least the mass of ten old-disc stellar 
particles) to the vertical velocity dispersion of all stars in 
a bin at (Ri) (a column in the plots) is represented. The 
stars migrating the most from an initial guiding radius tend 
to be colder in the z-direction than the average for all the 
stars at this initial guiding radius. 

We next investigated whether the migrators can be sub¬ 
ject to heating (or cooling) when they reach hotter (re¬ 
spectively colder) regions. In Fig. 18 we represent the dis¬ 
tribution of the ratio of the final vertical velocity disper¬ 
sion to the initial one for stars that have migrated by 
AR = (Rf) — {Ri) from their initial radius Ri in a time 
interval [ti,tf]. We observe that, especially in the inner re¬ 
gions of the disc, that is, inside the OLR, outward migra¬ 
tors tend to lower their velocity dispersion, while inward 
migrators tend to increase it. The effect increases with the 
amplitude of the variation (it is more visible for the extreme 
cases of migration). 

We now discuss how the final vertical velocity disper¬ 
sion of migrators compares to that of stars with the same 
final guiding radii (Rf). In Fig. we represent the distri¬ 
bution of the ratio of the velocit}Mispersion of migrators in 
a bin in (Rf) and (Rf) — (Ri) to the velocity dispersion of 
all stars of guiding radius (Rf). We observe different trends, 
depending, among others, on the location of the radius Rf 
under consideration. We start by discussing what occurs in¬ 
side the OLR at the final time. We observe that in the very 
inner regions, stars migrating inwards to reach a guiding ra¬ 
dius (Rf) have a higher velocity dispersion than the whole 
population of stars at (Rf). These stars must be taking 
part in the peanut-shaped bar with a high vertical veloc- 
ity dispersion . This agrees with the results reported by [DI] 


Matteo et al. (2014), who found that outside-in migrators 
participating in a boxy/peanut shaped structure tend to 
be kinematically hotter than in situ stars. At higher radii, 
however, the inward migrators tend to have lower velocity 
dispersions than the average at their final guiding radius. 
Outward migrators tend to have a slightly higher (~20%) 
vertical velocity dispersion than all stars at their final guid¬ 
ing radii, except for the most extreme migrators -those with 
the highest AR variations, found between the CR and the 
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Fig. 17. Ratio of vertical velocity dispersion of stars in bins in {R[) and {Rf) — {R[) to the vertical velocity dispersion of 
all stars in radial bin centred on (Ri). The shaded areas represent the bar ILR (blue), CR (red), and OLR (green) radii 
variation in the time-span of each plot. The red vertical lines are the average of the bar CR radius during the time-span 
of a plot. 



Fig. 18. Distribution of the ratio of final to initial vertical velocity dispersion of stars migrated by AR = (Rf) — {R[) 
from their initial guiding radius (Ri) in the time interval [ti, R]. The shaded areas represent the bar ILR (blue), CR (red), 
and OLR (green) radii variation in the time-span of each plot. The red vertical lines are the average of the bar CR radius 
during the time-span of a plot. 


OLR region. These extreme migrators can be significantly same final guiding radius. We note that this effect is only 
hotter (up to 50% more) than the whole population at the evident when a selection on the amplitude of the migration 
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Fig. 19. Distribution of the ratio of vertical velocity dispersion of stars in bins in (i?f) and (Rf) — {R[) to the vertical 
velocity dispersion of all stars with guiding radii centred on the final value (Rf). The shaded areas represent the bar ILR 
(blue), CR (red), and OLR (green) radii variation in the time-span of each plot. The red vertical lines are the average of 
the bar CR radius during the time-span of a plot. 



R[kpc] 

Fig. 20. Vertical velocity dispersion at t = 3 Gyr at R of all 
stars, non-migrators, or migrators, that migrated by more 
than 2 kpc in terms of radius or guiding radius from 1 to 
3 Gyr. 


AR is made. When the overall population of (inward or 
outward) migrators is analysed as a whole, the final verti¬ 
cal velocity dispersions of migrators and non-migrators are 
much more similar, which is because apart from the ex¬ 
treme cases of migration, many stars that have migrated 
by small AR tend to have a final kinematic similar to that 
of the overall population. We can see this in Fig.|^ In this 
figure, we have plotted the vertical velocity dispersion of all 
stars as a function of radius at t = 3 Gyr and the velocity 


dispersion of inward and outward migrators that are in a 
radius bin and that have migrated by more than 2 kpc in 
terms of radius or guiding radius since t = 1 Gyr, and the 
non-migrators defined as the stars of the radius bin that 
have not changed their radius or guiding radius by more 
than 2 kpc. The inward migrators have a slightly lower 
vertical velocity dispersion, the outward migrators have a 
slightly higher velocity dispersion, but the total velocity 
dispersion is very close to the one of the non-migrators, 
indicating a weak effect of the overall migration on the lo¬ 
cal velocity dispersion. We insist, however, on the fact that 
the most extreme migrators, in particular the most extreme 
outward migrators, can have dispersions significantly differ¬ 
ent (bO%}Jh;oni_thn averag e. While we agree with the results 
by Minchev et al. (2012a), therefore, that migration over¬ 
all contributes little to disc thickening, extreme migrators 
can depart from the average small increase of the heatin g 
found in Fig. 20 and found also by Minchev et al. (2012a), 
contributing vutn significantly higher vertical velocity dis- 
persions to the whole sample of migrating stars that end 
up at the same final guiding radius. 


6.2. Radial heating/cooling 

The radial velocity dispersion also increases on average with 
time, as can be seen in Fig.[^ We can follow the amplitude 
of the radial oscillations as well as the eccentricity of the 
trajectories of the stellar particles defined as 

-^max -^min /i-,\ 

^ ~ ~R ‘ 

-^max I JVjxiin 

The eccentricity varies between 0 for circular orbits 
and 1 for radial orbits. The mean eccentricity of the disc 
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Fig. 21. Time evolution of the radial velocity dispersion 
of the old stellar disc without the bulge (solid), the young 
stellar disc (dashed), and the gas disc (dotted-dashed). 




Fig. 22. Time evolution of the amplitude of radial oscilla¬ 
tions (top) and eccentricity (bottom) with time. 


stars increases with time. This increase is mainly caused by 
the growth of the bar: More and more stars are gradually 



Fig. 23. Change of eccentricity as a function of change of 
guiding radius from 1 to 3 Gyr (top) and 5 to 7 Gyr (bot¬ 
tom) . The colour scale represents the mass in each bin and 
is logarithmic. The solid black line is the average of the 
change in eccentricity by bin of change of guiding radius 
(i.e. the average of the corresponding column). 


trapped by the bar and gain eccentricity as their orbits be¬ 
come elongated. The evolution of the amplitude of radial 
oscillations (2 A(t) where A{t) is the semi-amplitude intro¬ 
duced in Sect. 4.2) and eccentricity e as a function of radius 
can be seen in Fig. The bar growth is visible from the 
increase of eccentricity and amplitude with time in the in¬ 
ner kpc. The amplitude of radial oscillations increases on 
average with time at all radii, which is consistent with the 
increase of radial velocity dispersion shown in Fig. and 
with the growing role of blurring with time, as discussed 
in Sect. 4.2 The eccentricity also slightly increases at all 
radii. 

It is interesting to investigate whether radial migration 
changes the radial amplitude and the eccentricity of the 
stars trajectories. In Fig. 23, we represent the change in 
eccentricity as a function of the change in guiding radius 
between t = 1 and 3 Gyr in the top panel and between t = 5 
and 7 Gyr in the bottom panel. The distribution is broad. 
The majority of stars only slightly change their eccentric¬ 
ity, but on average (black line), stars migrating outwards 
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(with positive (i?f) — (Ri)) tend to decrease their eccentric¬ 
ity while the opposite stands for stars migrating inwards. 
The inward effect increase of the eccentricity (especially 
visible in the bottom panel of Fig. 23 from 5 to 7 Gyr) is 
again dominated by the capture of stars by the central bar. 
The decrease of eccentricity does not necessarily mean that 
the orbits are circularised in terms of decrease in amplitude 
of radial oscillations because of the denominator in the ex¬ 
pression of the eccentricity, which can increase in the case 
of outwards migration. In Fig.j^we therefore represent the 
variation of the amplitude of radial oscillations and sepa¬ 
rate in the middle and right columns the disc into an ‘inner 
disc’ and an ‘outer disc’ by a cut between the CR and OLR 
bar radii at 10 kpc for the first time-interval and at 15 kpc 
for the second one, so that the effect of the bar growth on 
the stellar orbits is not present in the ‘outer disc’ plots. In 
the inner-disc plot of the time-span from 5 to 7 Gyr, a clear 
trend of an increase of the amplitude of the radial oscilla¬ 
tions is visible for inward migration, corresponding to stars 
migrating inwards and contributing to the bar. All the dis¬ 
tributions are again broad, with the largest variations of 
amplitude occurring for stars that do not migrate signifi¬ 
cantly. Stars that migrate outwards the most increase their 
radial oscillations amplitude on average, but the possible 
increase is limited compared with the rest of the distribu¬ 
tion. 

7. Conclusions 

We have studied stellar radial migration in a simulation of 
an Sb-type extended galactic disc. We confirmed the main 
role of corotations of the density resonances as seeds for 
radial migration in terms of both blurring and churning 
and observed a strong radial migration when there is a bar- 
sp iral resonance overlap; this is consistent with the results 
of Minchev & Famaey (2010); Minchev et al. (2011). 

We have quantified the effects of blurring and churning. 
Migration defined as the simple difference in galactocentric 
radii between two times can lead to significantly overesti¬ 
mating the fraction of migrators by churning - and also 
the spatial extension of this migration. The intrinsic na¬ 
ture of the epicyclic orbits and the increase in the ampli¬ 
tude of radial oscillations with time - both effects being 
included when migration is quantified by a change of in¬ 
stantaneous radius - can indeed give the impression of a 
much stronger migration than that really experienced by 
stars in the galaxy. This is particularly true, for example, 
for the fraction of migrators in the outer disc of our simu¬ 
lated galaxy. When blurring is excluded, the fraction of ex¬ 
treme migrators contributing to the outer disc population 
decreases by a factor between 2 and 8 (compare Fig. 13, 
top and bottom panels, for the time interval 1 to 9 Gyr). 

Whilst the spatial extent covered by migrating stars in¬ 
creases with time, our simulation suggests that migrators 
encounter barriers. In particular, stars migrating by churn¬ 
ing from corotation cannot cross the OLR, and vice versa, 
stars born beyond the OLR cannot reach the inner disc. The 
OLR region - defined as the region between the position of 
the OLR at the epoch of bar formation, and at the final 
epoch - is a transition region, the only region where some 
pollution between the inner and the outer disc is allowed. 

Even though our model is not intended to reproduce 
the Milky Way - the pattern speed of our simulated bar 
is lower than the pattern speed measured for the Galaxy 


(Gerhard 2011), and in consequence, the main resonances 
are located at much larger distances from the centre than 
those measured for the Milky Way, and the study assumes 
an initially already formed thin disc - we think that the 
previous result might help understand the puzzling nature 
of the outer Galactic disc and its significantly different stel¬ 
lar populations. It has recently been shown that the Milky 
Way outer disc followed a dif ferent chemical evolut i on his¬ 
tory than d id the inner disc ( Haywood et al.| 2Q13 Snaith 
et al.||2Q14 ). Regardless of the formation mechanisni of the 


outer regions of the Galaxy, stars there have been able to 
evolve independently of the inner disc. Our model suggests 
that the Galactic OLR and its location ( estimated to b e 
located slightly inside the solar radius, see Dehnen (2000|)) 
have a major role in explainin g this finding. The OLR in- 
deed acts as a barrier for gas (Gombes 1988), but - as we 
showed here - for stars as well, inhibiting migration from 
the inner to the outer disc and vice versa. This may explain 
why stars with i nner thin disc chemis try are not observed 
in the outer disc ( Haywood et al.|2Q13 ). This finding, if con¬ 
firmed with future studies, may also lead to doubts about 
the interpretation of the U-shape in age profiles, or inver¬ 
sion in colour-profiles found in the outer disc of external 
galaxies. We do not yet know whether our Galaxy has the 
U-profile in stellar ages that is sometimes observed in exter¬ 
nal galaxies like M33. When these inversions are observed 
outside the OLR position in barred galaxies, it is difficult 
to explain them in terms of strong migration from the in¬ 
ner disc, with the observed bar/spiral pattern speed, since 
they occur beyond their OLR. However, it is possible that 
previous bar/spiral waves have developed with lower pat¬ 
tern speeds, implying OLR and migration at larger radii. 
Another possible explanation is that the settlement of the 
Galactic disc in the outer regions permitted the formation 
of a significant nu mber of old stars in sit u, as proposed for 
the M ilky Way in Haywood et al. (2013) and Snaith et al. 

(iMa. 


Finally we have analysed the kinematics of migrating 
stars. 


We confirmed the results found by Vera-Giro et al. (2014) 
for spiral galaxies. Similarly to that case, also when migra¬ 
tion is mainly induced by a stellar bar, there is a provenance 
bias of migrators in terms of their kinetic state. The stars 
migrating the most from an initial guiding radius tend to 
be colder in the z-direction than the average of all the stars 
at the same initial guiding radiu^^ _ 


We also confirmed the results by Minchev et al. (2012a): mi¬ 
gration contributes li ttle to disc thickening, but we point 
out, consistently with Minchev et al. (2012b), that there is 
a trend of increasing vertical velocity dispersion with the 
extent of migration: the most extreme outward migrators 
that end up at a given final radius tend to also have the 
highest velocity dispersions when compared to the velocity 
dispersion of all the stars found at that final radius. Thus, 
while the overall effect of heating at a given radius is weak, 
we suggest that at a given radius, extreme outward migra¬ 
tors from the inner disc are identifiable as stars that have 
the highest velocity dispersions among those measured for 
stars of the same age, at the same radius. We recall that 
this signature is different from that induced by mergers, 
which could heat the outer disc enough for extreme migra¬ 
tors from the central pa rts of the galaxy to h ave a cooling 
effect on the outer disc ( Minchev et al.||2Q14 ). 


22 












































A. Halle et aL: Blurring, churning, and the outer regions of galaxy discs 



1 to 3 Gyr 


-10 -5 0 5 10 

<Rf >- <Ri > [kpc] 


< 




<Rf >-<Ri > [kpc] 



<Rf>-<Ri>[kpc] <Rf>-<Ri>[kpc] 


'.■J 


■ 

-■ iT 

Sj” 

Outer disc 


-10 -50 5 10 

<Rf >-<Ri > [kpc] 


Fig. 24. Change of radial oscillation amplitude as a function of change of guiding radius from 1 to 3 Gyr (top) and 5 to 
7 Gyr (bottom). The colour scale represents the mass in each bin and is logarithmic. The solid black line is the average 
of the change in amplitude by bin of change of guiding radius (i.e. the average of the corresponding column). 


Overall, our findings challenge the current view of the 
effect that radial migration from the inner disc may have in 
the outer regions of disc galaxies when only a main asym¬ 
metric pattern is present. This may fundamentally affect 
the understanding of stellar populations in bar-dominated 
galaxies, which we will investigate in future studies. 
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